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ABSTRACT 



The active galaxy Markarian 501 was discovered with air-Cerenkov telescopes 
at photon energies of 10 tera-electron volts. Such high energies may in- 
dicate that the 7 rays from Markarian 501 are due to the acceleration of 
protons rather than electrons. Furthermore, the observed absence of 7 ray 
attenuation due to electron-positron pair production in collisions with cosmic 
infrared photons implies a limit of 2 to 4 nanowatt per squaremeter per stera- 
dian for the energy flux of an extragalactic infrared radiation background at 
a wavelength of 25 micrometers. This limit provides important clues on the 
epoch of galaxy formation. 
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Gamma rays (7 rays) from cosmic sources impinging on Earth's atmo- 
sphere initiate electromagnetic showers in which the energy of the primary 7 
ray is imparted among secondary electron-positron pairs. The blue Cerenkov 
light emitted by the pairs in the atmosphere can be detected from the ground 
with optical telescopes triggering on the short (~ 1 ns) optical pulses. The 
technique has advanced considerably in recent years (1), and some surpris- 
ing discoveries have been made. Among them is the detection of the blazar 
Markarian 501 (Mrk 501) at energies above 10 TeV (1 TeV = 10 12 eV) (2). 

Blazars are remote but very powerful sources characterized by their vari- 
able polarized synchrotron emission. They are associated with radio jets 
(bipolar outflows) emerging from giant elliptical galaxies seen at small an- 
gles with the line of sight. Mrk 501 is ~ 3 x 10 8 light years from Earth but 
nevertheless produces a tera-electron volt 7 ray flux during outbursts that 
is many times stronger than that of the Crab nebula, a supernova remnant 
inside our Milky Way at a distance of only 6 x 10 3 light years. The radiation 
mechanism responsible for the 7 rays could be either inverse-Compton scat- 
tering of low-energy photons by accelerated electrons (3) or pion production 
by accelerated protons. In the latter case, the sources could be among the 
long-sought sources of cosmic rays, that is the isotropic flux of relativistic 
particles with differential number density (N) spectrum dN/dE oc E~ 2 - 7 (for 
energies E < 10 3 TeV) mainly consisting of protons and ions (4). 

Particle acceleration in astrophysics is typically observed to be associated 
with (collisionless) shock waves when a supersonic flow of magnetized ma- 
terial hits a surrounding medium. Examples of shock waves are shell-type 
supernova remnants (explosion of a massive star), plerions (pulsar wind), 7 
ray bursts (relativistic ejecta from the collapse of a compact stellar object), 
or the jets ejected from active galactic nuclei (collimated relativistic wind 
from the accretion disk around a supermassive black hole). 
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In the theoretical picture of shock acceleration, relativistic particles (pro- 
tons, ions, electrons) scatter elastically off turbulent fluctuations in the mag- 
netic field on both sides of the shock and thereby gain energy because of the 
convergence of the scattering centers (approaching walls). The acceleration 
time scale for the process can be written as t acc = £r g c/-u 2 where u denotes the 
velocity of the shock wave (c is the speed of light) and r g oc E/B denotes the 
radius of gyration of a particle with energy E in a magnetic field of strength 
B. The effects of shock obliquity, turbulence spectrum, and other unknowns 
are conveniently hidden in an empirical factor £ > 1. The most rapid (gyro- 
time scale) particle acceleration for relativistic shocks corresponds to £ = 1 

(5) . Balancing the acceleration time scale with the energy loss time scale 
due to synchrotron radiation t syn oc B~ 2 E~ X one obtains the maximum en- 
ergy of the electrons E max = 10 (f/10)-°- 5 (.B/3AiG)- - 5 (u/10 8 cm s" 1 ) TeV 

(6) . The observed 10 TeV 7 rays from the Crab nebula (7) and the observed 
synchrotron x-rays in shell- type supernova remnants (8) (corresponding to 
10 TeV electrons) require £ ~ 1 to 10. Because protons lose less energy, they 
can reach larger -EFax's than electrons and give rise to 7 ray emission even 
above ~ 10 TeV via pion production and subsequent pion decay. Although 
shock acceleration theory predicts that most of the cosmic rays are acceler- 
ated in supernova remnants (4), no definitive 7 ray signature has yet been 
discovered. 

It is commonly argued that the assumption of electron acceleration also 
suffices to explain the 7 rays from blazar jets such as Mrk 501 (9). Estimates 
of the magnetic field strength in the 7 ray emitting part of the jet in Mrk 501 
then yield values in the range B ~ 0.04 to 0.7 G. This magnetic field is much 
stronger than the one in supernova remnants and the associated stronger 
cooling of the relativistic electrons due to synchrotron energy losses reduces 
E max accordingly. The effect is almost compensated by the high shock wave 
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velocities in extragalactic radio sources which speed up the acceleration rate. 
Using radio interferometry, shock wave velocities close to the speed of light 
have been inferred corresponding to typical bulk Lorentz factors in the range 
Tjet = (1 — /3 2 ) -0 - 5 ~ 2 to 10 (f3 — u/c) with a few cases of still higher values 
(10). Due to the alignment of the jet axis and the line-of-sight in Mrk 501 
superluminal motion has not been observed. With u — c, £ — 10, and tak- 
ing into account equal synchrotron and inverse-Compton losses, one obtains 
-Smax ~ 4 rj Ct (-B/G)~ ' 5 TeV from the balance between acceleration gains 
and energy losses. The additional factor rj ct takes into account the boost 
in energy due to the relativistic bulk motion. Therefore, electron maximum 
energies of ~ 10 TeV as required for Mrk 501 (at least 5-8 TeV are required 
for the similar blazar Mrk 421 (11)) are formally allowed, but one certainly 
has to push the theory to its limits and this raises a number of concerns: (i) 
the multi-TeV spectrum should show considerable curvature due to the (so- 
called Klein-Nishina) decrease of the scattering cross section when the energy 
of the scattered photon approaches E max and due to the onset of electron- 
positron pair production (the observed multi-TeV spectrum is consistent with 
a smooth power law), (ii) the ratio between the 7 ray and synchrotron (si- 
multaneous) luminosities depends sensitively on the jet Lorentz factor Fjet 
and therefore requires fine-tuning (both nearest bright blazars, Mrk 421 and 
Mrk 501, show a similar 7-to-x-ray luminosity ratio), (iii) the magnetic field 
pressure turns out to be much lower than the relativistic electron pressure 
in the electron acceleration models which seems inconsistent with the shock 
acceleration mechanism (the turbulent magnetic field is responsible for push- 
ing the electrons back and forth across the shock), larger values of B are also 
expected from the adiabatic expansion of a magnetically collimated jet (the 
observed S-field at the tips of the jet in Cygnus A is consistent with adiabatic 
expansion (12)), and (iv) larger values of rj et could ameliorate the problem 
that E max ~ 10 TeV, however, one would run into a serious problem with 
unification models of active galaxies if Fjet > 10 would be the rule rather 



Karl Mannheim 6 



than the exception (13) (number of required host galaxies would exceed the 
number of known radio galaxies). 

A natural solution to the problem is to assume that the 10 TeV 7 rays 
are due to pion production from accelerated protons. The balance equa- 
tion between energy gains and losses for protons yields maximum energies 
of ~ 10 6 TeV and short variability time scales t var > 10 5 (^/r jct )( J B/G)- 1 s 
in Mrk 501 (14)- The relativistic proton energy loss is dominated by photo- 
production of pions in collisions with low-energy synchrotron photons (origi- 
nating from accelerated electrons). Collisions of the accelerated protons with 
matter are negligible due to the low matter density in relativistic jets (unless 
a high-density target moves across the jet (15)). The 7 rays from the decay 
of the neutral pion (far above the observed range of energies) are subject to 
pair creation in further collisions with the low-energy synchrotron photons 
(7 + 7 — > e + + e~). This initiates an electromagnetic cascade shifting the 
average photon energy to the TeV range and below. A model based on the 
combined acceleration of protons (7 rays) and electrons (radio-to-x-rays) - 
coined the proton blazar model (16) - was fitted to published data of Mrk 501 
in order to obtain a prediction of its multi-TeV spectrum (17). Data from 
1995 and earlier was available for the analysis and covered the radio-to-x-ray 
wavelength range including a flux limit above 100 MeV and an integral flux 
above 300 GeV (for details see references in (17)). The published flux values 
showed considerable variability in the optical-to-x-ray range and the model 
spectrum was therefore fitted to match the time-averaged spectrum (from 
the fit one obtains B = 37 G and r je t = 10). The predicted multi-TeV spec- 
trum is shown in Fig.l and compared with the data obtained from recent 
(1996, 1997) air-Cerenkov observations. The fairly robust spectral slope of 
the model spectrum fits the observations very well, whereas the absolute flux 
normalization is somewhat too low. Considering that the sub- TeV (350 GeV) 
flux has been reported to increase from 1995 to 1997 (9), the agreement is 
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actually rather impressive if one scales up the model spectrum accordingly. 
Contemporaneous multi-wavelength observations of blazars such as Mrk 501 
will be important to discriminate between electron-based and proton-based 
models for the 7 ray emission from them. Generally, electron-based models 
require larger values of rj ct and lower values of B than proton-based models 
to obtain high-energy 7 rays. 

There is a further hint that proton acceleration might be important. Un- 
resolved blazars are the most probable source population to produce the 
observed extragalactic 7 ray background between 10 MeV and 10 GeV (18). 
The energy density of this 7 ray background is ~ 4 x 10" 6 eV cm -3 . A sim- 
ilar value is found for an extragalactic flux of protons (with dN/dE oc E~ 2 
differential spectrum between 10 9 eV and 10 20 eV) providing all of the ob- 
served cosmic rays with energies above 10 18 - 5 eV (the so-called "ankle" above 
which the very high energy differential cosmic ray spectrum oc E~ 3 flattens) 
(19). On the assumption that the 7 rays are from proton (p) acceleration, 
the comparable energy densities result from simple decay kinematics: Photo- 
production of neutral pions (n°) (p + 7 — > n° +p) and their subsequent decay 
gives rise to 7 rays (subject to electromagnetic cascading) which carry ~ 1/5 
of the proton energy. Charged pions (p + 7 — »• n + + n) are produced with 
approximately the same rate and give rise to neutrons (n) which carry ~ 4/5 
of the accelerated proton energy. Since the neutrons do not scatter off the 
magnetic field fluctuations responsible for the acceleration and storage of the 
charged particles in the blazar jet, they must escape the accelerator (energy 
losses are small and modify the neutron spectrum only at the very highest 
energies). Time dilation allows the most energetic neutrons to leave the host 
galaxy freely before /9-decay occurs (for protons there would be adiabatic 
losses due to the magnetic field in the host galaxy). Hence the neutron (and 
after /5-decay proton) luminosity is equal to the 7 ray luminosity within fac- 
tors of order unity. An extragalactic origin of the highest energy cosmic rays 
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is indeed suggested by the absence of an enhancement of the cosmic ray flux 
toward the Galactic disk (20) and by the change in chemical composition 
from heavy (protons and ions) to light (protons) above 10 18 5 eV (21). The 
ultimate challenge to the hypothesis is the measurement of the high-energy 
neutrino (z/) flux associated with the charged pion decay (n^ — > e ± + 3u). 
The energy density in these multi-TeV neutrinos would be of the same order 
of magnitude as that in extragalactic cosmic rays and 7 rays. Their mea- 
surement therefore constitutes an experimentum crucis within reach for the 
planned cubic kilometer underwater(ice) detectors (22). 

Although 7 rays are known from laboratory experiments for their pene- 
trating power, propagation over intergalactic distances is not without hurdles. 
A diffuse isotropic infrared background (DIRB) was produced when the first 
galaxies formed. Massive stars in early galaxies produced large amounts of 
dust in their winds reprocessing the visual and ultraviolet light from the 
stars into infrared light. By colliding with these ample infrared photons, 
7 ray photons can disappear and turn into electron-positron pairs (23-25). 
The most numerous infrared photons above threshold for pair production 
with 10 TeV 7 rays have wavelengths ~ 25 //m. The mean free path (A 77 ) 
for pair creation at multi-TeV energies is of the order of the distance (d) of 
Mrk 501. The exact value depends on the DIRB which is difficult to mea- 
sure directly owing to the presence of zodiacal light and galactic cirrus clouds. 

One can use the observed power law spectrum (2) to put a limit on the 
maximum allowed pair attenuation assuming that the observed power law 
is the unattenuated spectrum emitted by the source (consistent with the 
proton-based model). In general, only very contrived intrinsic spectra would 
look like a smooth power law after the quasi-exponential attenuation. The 
maximum allowed deviation from the power law (1 — exp[— d/A 77 ]) is taken 
to be the size of the statistical error bar at 10 TeV yielding an optical depth 
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t 77 = d/A 77 < 0.7. This limit can be relaxed by a factor not larger than ~ 2 
admitting for weakly absorbed spectra which still approximate a power law 
(see the dashed line in Fig.l). There is some dependence of the attenuation 
on the shape of the DIRB spectrum. Useful models for the spectral shape can 
be found in (23-25) and yield a similar limit for the 25 /zm DIRB normaliza- 
tion ul u (25 am) < 2 to 4 nW m -2 sr^ 1 . The absence of 7 ray attenuation in 
Mrk 501 is consistent with no contribution to the DIRB other than from the 
optically selected galaxies for which one expects ~ 10% of their optical emis- 
sion to be reprocessed by warm dust yielding ul u (25fim) ~ 1 nW m -2 sr _1 
(26), but would also allow a DIRB stronger by a factor 2 to 4. A DIRB 
of at least ~ 3 nW m -2 sr _1 is suggested by faint infrared galaxy counts 
and indicates contributions from dust-enshrouded galaxies at redshifts of 
z ~ 3 — 4 (24)- Electron-based models for the 7 ray emission from Mrk 501 
(9) predict deviations from a power law in the multi-TeV range even with- 
out external attenuation and therefore impose an upper limit on the DIRB 
below the lower limit from faint infrared galaxy counts. If both methods to 
estimate the DIRB (deviations from a power law spectrum in the multi-TeV 
range, faint infrared galaxy counts) use correct assumptions, a cutoff in the 
7 ray spectrum of Mrk 501 must be present in the energy range 10 to 30 TeV. 
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Figure captions 

Fig.l: Differential TeV spectrum of Mrk 501. Thin solid line: proton blazar 
model prediction based on archival data from 1995 and earlier (17) taking 
into account intergalactic attenuation adopting a DIRB spectrum based on a 
cold+hot dark matter galaxy formation model from (23) and normalized to 
vl v (2hnm) = 1.0 nW m -2 sr" 1 (a Hubble constant of H Q = 75 km s" 1 Mpc" 1 
is assumed throughout the paper). Thin dashed line: same without inter- 
galactic attenuation. Solid thick lines: model spectrum scaled up by factors 
2.5 and 20 corresponding to the increase in the mean sub- TeV flux from 1995 
to 1996 and spring of 1997, respectively (9). Open circles: HEGRA (High En- 
ergy Gamma-Ray Astronomy) CT1 observation March-August 1996 (220h) 
(27). Solid circles: HEGRA IACT observation March- April 1997 (26. 7h) (2) 
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